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ABSTRACT 

In this paper, four sets of evolutionary models are computed with different values of the mixing 
length parameter a p and the overshooting parameter S ov . The properties of the convective cores 
and the convective envelopes are studied in the massive stars. We get three conclusions: First, 
the larger a p leads to enhancing the convective mixing, removing the chemical gradient, and 
increasing the convective heat transfer efficiency. Second, core potential <f> c = M c /R c describes 
sufficiently the evolution of a star, whether it is a red or blue supergiant at central helium ignition. 
Third, the discontinuity of hydrogen profile above the hydrogen burning shell seriously affect the 
occurrence of blue loops in the Hertzsprung Russell diagram. 

Subject headings: stars: evolution — stars: interiors—stars: supergiants 


1. INTRODUCTION 


Massive stars have convective cores and ra¬ 
diative envelopes. Convection can mix different 
chemical elements in the convection zone, and 
bring heavy elements generated by nuclear re¬ 
action region to the stellar surface. It is also 
an important energy transport mechanism. Lo- 
cal mixing length th eory (MLT) developed by 
Bohm-Vitense ( 19581) is the most widely used 
description for convection in stellar model. The 
mixing length is proportional to the pressure scale 
heig ht l = a„F n , wher e a n is a free parame - 
ter. IPasetto et al. ( 20141) and Magic et al. ( 2014 ) 
showed that the mixing length parameter a p was 
of great importance in determining the convective 
energy transport, and directly affected the stellar 
interior structure of temperature. However, from 
the theoretical point of view the knowledge of a p 
is still a long standing problem. 

Convective overshooting belongs to convection 
and also affects energy transfer and mass mixture 
in stars. Classical treatment of the convective 
overshooting i s based o n the non-local mixing 
length theory (, Zahn 19911 ). The predicted proper- 
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ties of stellar evolution are highly sensitive to the 


sequence and core helium burning phases (Maeder 

1975, 

19761 Stothers & Chin 1981: Bertelli et al. 

lO 

00 

T—1 

Noels etalj 20ld Chen & Han 20oi 2003| 

Montalban et al. 20131). Moreover. Bressan et al. 

(2014) and Ding & Li 

2014) showed that con- 


vective overshooting and semiconvection greatly 
affected the extent and the rate of mixing in the 
outside deeply convective envelopes. 

Red supergiants (RSGs) are evolved, massive 
(10 — 30 i) co r e helium burning star s (e.g. , 


Mevnet fe Maedcil l2000t Eldridge et al. 20081: 


Brott et al. 2011 : iDavies et al. 20131 ). It 


repre¬ 
sents a key phase in the evolution of massive stars. 
However, convection overshooting in the evolution 
of stellar models from blue to red in the post main 
sequence phase and from red to blue in the blue 
loop phase have not been studied systematically. 


In this paper, we attempt to investigate the 
properties of the convective cores and the con¬ 
vective envelopes of massive stars under different 
convection and convective overshooting assump¬ 
tions. Special attention is given to the convective 
overshooting at the bottom of convective envelope 
during the RSG phase. §2 describes the evolution 
code, and gives the input parameters for the 
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Sov 

Case 1 

1.5 

0.0 

Case 2 

1.5 

0.12 

Case 3 

2.0 

0.0 

Case 4 

2.0 

0.12 


Table 1: Four sets of evolutionary models. 


modeling. §3 gives the main results and discussion 
in detail. In §4 we summarizes our conclusions. 

2. MODEL 

We employ an updated version of stellar evo¬ 
lution code ev, which is based on the evolution 


program of Eggleton ( 197ll 1972 . 1973 ). In this 
paper, we use 499 meshpoints in order to be more 
accurate. Evolutionary models are computed for 
10, 15, 20, 30, 40, 60 Mq stars. For each mass, 
four sets of evolutionary models from the main 
sequence (MS) to the end of core helium burning 
phase are computed (see Table 1). The first and 
third sets are computed without a certain amount 
of overshooting and computed with radios a p = 
1.5 and 2.0, respectively. The second and fourth 
sets are co mputed with convec t ive overshooting 


S ov = 0.12 ( Schroder et al.l 1997; P ols et al.l ll998) 


and computed with radios a p = 1.5 and 2.0, 
respectively. We neglect the effects of rotation, 
mass loss and magnetic fields. 


3. RESULT 
3.1. H-R diagram 

Figure 1 shows four cases of evolutionary tracks 
in the Hertzsprung- Russell diagram (HRD) for 
each mass. We can see that t he stars evolve t oward 
lower effective temperature. lEldridge et al.l ( 20081) 
took logT e ff < 3.66 as a condition to be a RSG 
during the central helium burning stage. First we 
compare cases 1 and 2 for each mass. It can be 
noted that convective core overshooting produces 
an extended track during MS phasejmd case 2 has 
a larger luminosity than case 1. iMaeder &; Mevnet 
( 19871) suggested that models with an MS ter¬ 
mination in the range of spectral type 50 or 
later (logT e g < 4.47), the adopted convective 
overshooting model produce an MS widening. On 
one hand, convective core overshooting makes the 
nuclear reaction region larger, increases the MS 



Fig. 1. — Evolutionary tracks of 10, 15, 20, 30, 
40, and 60 Mq stars on the H-R diagram, from 
the zero-age main sequence to the exhaustion of 
central helium. 


lifetimes and makes the chemical discontinuity 
displace farther from the stellar core. On the other 
hand, the extension of mixing increases more fresh 
hydrogen in the convective core and enhances the 
energy production. The obtained conclusions are 
also applicable to cases 3 and 4. 

Second, we check cases 2 and 4 for each mass. 
We find that case 4 has a higher effective tem¬ 
perature than case 2 during the central helium 
burning phase. It indicates that conducting heat 
of convection is effected deeply by a p . In the view 
of fluid dynamics, a larger value of a p means that 
buoyancy can do more work on the convective 
cells. Thus, the degree of convective mixing is 
enhanced, and the efficiency of convective heat 
transfer is increased. 

3.2. Internal structure 

3.2.1. Lifetimes 

Four cases of the internal structures of 20 
and 30 Mq stars are showed in Figure 2. It is 
evident that core overshooting significantly puts 
off the occurrence of the helium burning phase 
and increases the MS lifetimes. iMaeder &; Mevnet 
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Fig. 2.- The internal structure of 20 and 30 M 0 stars, which are from the MS to the exhaustion of central 
helium. Filled regions represent the fully convective zones. 


(1987) found that lifetime depended on the ratio of 
q cc M/L, where q cc is the initial core mass fraction, 
M is stellar mass and L is stellar luminosity. In 
cases 2 and 4, there are more massive convective 
cores at the MS phase, and the thermonuclear 
reactions increase the MS lifetime tn- However, 
during helium burning phase the lifetime tne is 
decreased by convective core overshooting. Due 
to the evolutionary tracks of models still stay at 
a higher luminosity. As a result, convective core 
overshooting reduces the tne/tn ratio. 


3.2.2. Convective envelope 


Our models show that during RSG phases the 
stars possess deeply convective envelopes, which 
locate in the superadiabatic layer at the top rather 
than the hydrogen burning shell at the bottom. 
Lai fe Lil ( 2012 ) indicated that diffusive mixing in 


the overshooting region would lead to the increase 
of opacity and radiative temperature gradient 
V ra d- As a result, a convective stable region 
below the base of the convective envelope will be 
converted to a convective unstable region. 


It should be noticed that case 1 has a larger 
convective envelope than case 3, and case 2 has 
a larger convective envelope than case 4. This 
indicates that the larger a p produces more efficient 


mixing in the chemical gradient region, leads 
to a negligible chem ical gradient and suppr esses 
the semiconvection Stothers fc Chinl (1985) and 
Bressan et al.l ( 198ll) found that core overshooting 


reduced convective instability in the envelope. 


It is striking to notice that both cases 1 and 
3 have a small intermediate convective zone, just 
above the convective core, in the early stage of 
core helium burning. The reason is that the 
model without core overshooting leads to a higher 
hydrogen abundance outside the convective core. 
As a result, opacity and V ra d will be larger. So a 
semiconvection as well as full convection observed 
in the early core helium burning phase. 


3.2.3. Stellar radius 

Figure 3 shows the behaviour of radius R (in 
solar units) as a function of time. We can see that 
case 1 has a larger radius than case 3, and case 
2 has a larger radius than case 4 during helium 
burning phase. The radius of a star is defined as 



Enhancing the efficiency of convection (a p ) will 
increase the effective temperature of a star. Hence 
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Fig. 3.- The variations with time of the radius. 


radius decreases at a given luminosity. lLai fc Li 
(120121 ) also suggested that the expansion of con¬ 
vection in the stellar envelope increase convective 
heat transfer efficiency, and lead to a smaller 
stellar radius and a higher effective temperature. 


3.2.4■ Core mass 

In figure 4, for each mass the evolution of the 
core mass as a function of the time is presented. It 
should be noticed that cases 1 and 3 overlap, and 
cases 2 and 4 overlap. It shows that the different 
value of a p have no effect on the mass of core. 


3.3. Core potential 


When does a ma ssive star enter into RSG 
phase? Lauterborn et al. (Il97ll) suggested that 
the core potential <j) c = M c /R c describe suffi¬ 
ciently the evolution of stellar models from blue 
to red in the post main sequence phase. M c and 
R c are the mass and radius (in solar units) of the 
core, which extend from the centre up to the base 
of the hydrogen burning shell. 


We follow the work of iMowlavi & Forestini 
( 1994 ). When a star closes to the Hayaslri line, a 
critical potential cj> i cc (defined by log(T e ff) = 3.6) 
can be found. When </> c exceeds the critical value, 
the star becomes a red (super)giant. Furthermore, 


Fig. 4.- The variations with time of the mass of 
central convective hydrogen and helium core. 


according to the work of lEldridge et al.1 ( 200 8!), 
we define the critical potential face (defined by 
log(T’eff) = 3.66) above which the star enter into 
RSG phase. Figure 5 shows the core potential 
changes with the time in 20 and 30 M© stars. 


We can see that <f>i cc and <j> 2 CC show a good 
agreement. In cases 1 and 3, the </> c remain a 
long time under the <p cc . The star is still in the 
blue when central helium ignition. The reason is 
that the hydrogen burning shell grows rapidly in a 
region of smoothly varying hydrogen abundances. 
When <j> c reaches 0 CC , the star becomes a red 
supergiant. In cases 2 and 4, <p c rapidly exceed 
(f> cc . This result can be understand by considering 
that the convective overshooting makes a larger 
core during the MS phase. The star begins central 
helium burning as a red supergiant. 


3.3.1. Blue loop 


However, when the hydrogen burning shell 
encounters a discontinuity of hydrogen profile 
during the central helium burnin g stage, the issue 
will r elate t o the blue l oop (jStothers fc Chin 
1981 . 1 199 it lAlongi et al.1 1991 ). Our models 


develop huge convection envelopes during RSG 
phases. Downward convective overshooting at 
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Fig. 5.- The variations with time of the core potential in 20 and 30 Mq stars. The dotted line represents 
the </> ice- The vertical full lines refer to the See¬ 


the bottom of convective envelope will drive the 
hydrogen discontinuity deeper into the interior. 
When the hydrogen burning shell moves outward 
and encounters the discontinuity, fresh hydrogen 
will supply more nuclear fuels. So the energy 
production is enhanced, and the luminosity as well 
as the effective temperature is inc reased. Th e star 
then becomes a blue su pergia nt. [Lai fc Lil J2012h 
and Lauterborn et al. ( 197 lh showed that the 
occurrence of blue loop critically depended on the 
hydrogen profile just above the hydrogen bu rning 
shell. However, Huang fe Weigertl ( 1983 1 and 
iBertelli et al.1 ( 1985tl found that convective core 
overshooting during MS phase could significantly 
suppres s the sub s equen t development of the blue 
loop, [x u fe Lil (2004h showed that convective 
core overshooting took away the chemical gradient 
profile outside the stellar core, as well as the 
abundance discontinuity formed later by down¬ 
ward penetration of the convection envelope. 


4. Conclusions 

In this paper, we investigate the properties 
of convective cores and convective envelopes of 
massive stars. Different values of the overshooting 


parameter 6 and the mixing length parameter a p 
result in the different theoretical predictions. The 
main conclusions are summarized as follows: 

1) Convective overshooting makes a larger con¬ 
vective core, increases the MS lifetime and extends 
MS band in the HRD. 

2) Increasing the mixing length parameter a p 
leads to enhancing the convective mixing in the 
convective region and increasing the convective 
heat transfer efficiency. Thus, the star has a 
smaller chemical gradient, a higher effective tem¬ 
perature and a smaller stellar radius. 

3) The core potential </> c = M c /R c describes 
sufficiently the evolution of a star when it moves 
from blue to red. Without convective over¬ 
shooting, cf > c slowly reaches <fi cc and the star 
becomes a red (super)giant. The reason is that 
the hydrogen burning shell grows in a region of 
smoothly varying hydrogen abundances. With 
convective overshooting, </> c rapidly exceed <j) cc and 
then the star becomes a red (super)giant. Because 
a larger core is produced during the MS phase. 

4) The discontinuity of hydrogen profile above 
the hydrogen burning shell is play an important 
role in the occurrence of blue loops in the HRD. 
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